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SECTION  1.  SUMMARY 

1 0 BACKGROUND 

The  U.S.  Arniy  Yuma  Proving  Ground  (YPG)  air  delivery  test  mission 
includes  both  air  drop  and  air  transportability  tests  under  the  general 
connotation  of  Air  Delivery  Tests  are  conducted  using  three  separate 
drop  zones  for  dynamic  testing  from  aircraft,  and  a controlled  impact 
facility  for  static  drop  testing.  The  drop  zones  are  designed  for 
explosive  devices,  nonexplosive  devices,  and  personnel. 

The  Naval  Air  Facility  at  El  Centro,  California  (93  miles  from  YPG) 
offers  facilities  to  supplement  and  extend  the  YPG  capabilities.  YPG 
also  conducts  tests  in  nearby  Senator's  Wash  Reserve’ r on  the  Colorado 
River  for  the  conduct  of  water  impact  tests. 

In  addition,  YPG  provides  both  personnel  and  instrumentation  to 
support  air  delivery  tests  at  other  installations  such  as  Fort  Bragg, 
North  Carolina;  Fort  Carson,  Colorado;  and  White  Sands  Missile  Range, 

New  Mexico. 

The  YPG  air  transportability  ( internal /external ) mission  is  becoming 
increasingly  important  and  its  workload  increasing  accordingly.  Major 
tests  in  this  area  include  the  FAMECE  (Family  of  Military  Engineers 
Contractor  Equipment)  Test  Series,  the  UET  (Universal  Engineer  Tractor) 
Test  Program,  the  Battery  Control  System  Test  Program  aiyl  the  Armor 
Personnel  Carrier  Test  Series. 

The  existing  vibration  data  analysis  is  done  by  two  different 
groups  at  YPG.  The  accelerometer  instrumentation  and  playback  is  han41e4 
by  the  Electronics  Section  with  the  use  of  the  usual  airborne  and 
laboratory  magnetic  tape  recorders,  oscillographs  for  time  history 
editing,  and  a narrow-band  spectral  analyzer.  The  detailed  data  analysis 
is  accomplished  by  the  Data  Reduction  Section  on  an  as-available  class 
service,  which  usually  takes  several  weeks. 

Even  longer  times  are  involved  when  air  transportability  tests 
terminate  at  bases  other  than  YPG,  and  the  data  tapes  must  be  brought 
back  to  YPG  before  any  data  reduction  analysis  can  be  performed. 

The  long  analyzer  computation  time,  the  resulting  lack  of  data 
analysis  flexibility,  and  the  desire  to  give  the  project  engineer  a 
choice  in  data  formats  has  led  to  the  development  for  a small  dedicated 
system  that  could  be  operated  in  the  field. 

A Methodology  Investigation  of  an  airborne  recording  system  devel- 
oped by  the  Air  Delivery  Engineering  Branch,  YPG,  was  conducted  from 
November  1976  through  December  1977. 


2.0  OBJECTIVE 


This  report  describes  the  hardware  and  software  implementation  of 
an  airborne  recording  system  for  recording  and  processing  accelerometer 
vibration  spectral  data  from  air  delivery  and  transportability  tests  of 
AnT\y  materiel.  The  system  is  based  on  a programmable  calculator  and 
record/reproduce  instrumentation.  An  HP9825A  calculator,  HP9871A 
impact  line  printer,  and  interface  electronics  were  integrated  with  an 
analog  14-channel  magnetic  tape  recorder/ reproducer  to  digitally  process 
and  print/plot  acceleration,  power  spectral  density  and  displacement 
verses  frequency  from  several  tri axial  accelerometers.  The  system  is  an 
outgrowth  of  a long  standing  requirement  of  the  Air  Delivery  Engineering 
Branch  of  YPG  for  a portable,  stand-alone  system  for  air  transportability 
testing  of  am\y  materiel  test  items  to  AR  70-41  standards. 

3.0  SYSTEM  DESCRIPTION 

3.1  Summary 

The  programmable  calculator-based  vibrational  analysis  system  can 
best  be  understood  by  studying  the  block  diagram  illustrated  in  Figures 
1 and  2. 

First  there  is  a data  acquisition  phase  (Fig.  1)  where  the  outputs 
from  the  tri axial  accelerometers  are  recorded  on  magnetic  tape  for  later 
data  analysis.  At  the  same  time  IRIG-B  time  code  and  any  voice  annuncia- 
tion signals  are  also  recorded  for  the  purpose  of  knowing  where  to  index 
into  the  tape  to  find  the  data  of  interest.  All  of  this  equipment  is 
capable  of  operating  off  28-volt  dc  aircraft  power. 

Next,  there  is  a data  analysis  phase  (Fig.  2)  where  the  analog  data 
from  the  magnetic  tape  reproducer  is  digitized  for  input  into  the  program' 
mable  calculator.  Furthermore,  the  system  makes  use  of  a tape  search 
subsystem  to  find  starting  and  ending  points  of  data  channels  for  an 
operator's  predetermined  time.  This  subsystem  uses  the  IRIG-B  time  code 
channel  recorded  in  the  data  acquisition  phase  for  this  information. 

The  digitized  accelerometer  data  read  into  the  calculator  is  then 
transformed  into  the  frequency  domain  by  use  of  a Fast  Fourier  transform 
program.  This  program  outputs  power  spectrum  versus  frequency  informa- 
tion to  several  plotting  routines  for  different  plotting  formats.  These 
plots  are  then  outputted  to  an  impact  line  printer  for  display  to  the 
project  engineer. 
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3.2  Hardware 


3.2.1  Calculator 


The  Hewlett-Packard  9825  calculator  and  its  peripheral  as  configured 
for  this  application  are  shown  in  the  Figure  3.  The  calculator  is 
programmed  in  HPL  which  is  a Hewlett-Packard  version  of  multiple  statement 
per  line  BASIC  language.  HPL  executes  faster  than  standard/basic.  The 
calculator's  work  in  conjunction  with  24K  words  of  16-bit  read/write 
memory  and  a number  of  read  only  memory  (ROM)  units.  The  ROMs  required 
for  this  application  are  Advanced  Programming  ROMs.  They  are  used  as 
follows:  FOR  NEXT  statements.  String  for  character  string  manipulation. 
Matrix  for  array  manipulation.  General  I/O  and  Extended  I/O  for  buffered 
I/O,  and  HP-IB  control.  The  calculator  includes  an  alphanumeric  keyboard 
and  LED  display  for  convenient  operator  Interaction. 


FIGURE  3.  Configuration  of  the  Hewlett-Packard  9825.  calculator  and 
peripheral  equipment. 
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3.2.2  Printer 


An  HP9871A  impact  printer  is  the  peripheral  device  utilized  for 
displaying  the  computed  spectral  functions.  This  printer  is  a character- 
serial,  impact  printer  capable  of  performing  all  normal  typewriter  and 
printer  functions,  plus  plotting  functions.  All  printer  operations  are 
controlled  from  the  calculator  through  an  HP-IB  interface  bus.  This  bus 
is  the  same  as  the  IEEE  standard  488  general  instrumentation  bus  and  is 
used  to  provide  the  system  the  capability  for  future  expansion.  All 
plotting  commands  and  subroutines  needed  for  this  printer  are  included 
in  a tape  cartridge  of  general  utility  routines. 

3.2.3  Analog-to-Diqital  Converter 

A custom  made  analog-to-digital  conversion  module  that  takes  the 
analog  accelerometer  data  outputs  from  the  magnetic  tape  reproduction 
unit  was  built  into  the  system.  This  module  changes  the  analog  voltage 
to  a parallel  digital  signal  that  the  programmable  calculator  can  read 
and  operate  on  through  its  16-bit  duplex  interface.  This  circuit  consists 
of  a high-speed  sample-hold  amplifier  which  is  designed  to  acquire  and 
hold  up  to  a +10-volt  analog  signal  to  an  accuracy  of  jj3.01%  of  full 
scale  range  in: 

a.  5.0  microseconds  for  a 20-volt  step 

b.  4.5  microseconds  for  a 10-volt  step 

The  analog-to-digital  converter  (Fig.  4)  uses  successive  approximations 
with  a 12-bit  resolution  to  output  parallel  digital  data  which  is  TTL 
compatible.  This  means  that  upon  receipt  of  an  external  CONVERT  command, 

the  previous  data  sample  is  cleared  from  the  output  register  and  each 

bit  is  then  successively  compared  against  the  amplitude  of  the  input 
signal  and  is  either  held  as  a "0“  or  turned  on  as  a "1"  until  all  bits 
have  been  tried.  The  parallel  data  output  is  not  available  for  transfer 
to  the  calculator's  interface  until  the  STATUS  output  changes  from  logic 
"1"  to  logic  "0".  The  external  CONVERT  command  is  generated  form  a 
100  kHz  crystal  oscillator  and  CMOS  frequency  divider  chain  illustrated 
in  Figure  5.  The  final  2.5  kHz  signal  from  this  chain  is  fed  to  a mono- 
stable circuit  which  narrows  the  pulse  width  for  the  CONVERT  command. 

The  CONVERT  command  is  gated  with  the  valid  interval  signal  from  the 
tape  search  system  so  that  the  A to  D converter  outputs  data  only  between 
the  preset  START  and  STOP  times. 

The  STATUS  output  from  the  A/D  converter  is  connected  to  the  peri- 
pheral flag  line  on  the  16-bit  duplex  interface.  This  "tells"  the 
interface  to  read  the  converted  digital  value  from  the  A/D  output  reg- 
ister to  the  interface's  input  register.  The  STATUS  output  from  the  A/D 
converter  is  also  connected  to  the  peripheral  status  line  on  the  16-bit 
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A/D  CONVERTER 


FIGURE  4.  Sample-hold  amplifier  with  analog-to-dlgital  converter  circuit 
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duplex  interface.  This  is  done  so  that  the  calculator  can  check  this 
line  under  software  control  and  read  the  value  from  interface's  input 
register  into  its  own  memory  buffer. 

3.2.4  Tape  Search  Subsystem 

A Tape  Search  and  Control  Unit  is  designed  into  the  system  to 
provide  the  necessary  controls  to  automatically  search  the  data  tape  to 
locate  preselected  START  and  STOP  times.  This  device  is  used  with  the 
time  code  reader  to  accomplish  this  task.  This  START  and  STOP  time  is 
used  later  to  control  the  length  and  location  from  the  tape  of  the 
digitized  accelerometer  data.  Panel  lights  provide  front  panel  indica- 
tion of  the  status  of  the  tape  transport  as  well  as  the  entire  system. 
Manual  operation  of  the  tape  transport  is  also  possible  with  the  pushbut- 
tons. Adjustments  are  provided  to  optimize  tape  travel  during  high- 
speed search.  For  example,  the  search  operation  is  implemented  based  on 
time  frame  detection  and  on  real-time  durations.  Since  recording  speeds 
vary  with  respect  to  search  speeds,  tape  distances  traveled  also  vary 
drastically. 

Tape  transports  starting  and  stopping  characteristics  also  vary. 
During  search,  it  is  necessary  that  the  tape  be  stopped  ahead  of  start 
time  a sufficient  distance  so  that  detection  of  start  can  be  made  at  a 
controlled  speed,  always  in  the  same  direction.  It  is  undesirable  to 
begin  playback  too  far  in  front  of  START  time  since  valuable  time  will 
be  wasted  before  reaching  START  time.  Likewise,  if  playback  is  initiated 
too  close  to  START  time,  it  could  be  missed  due  to  transport  start  up 
times  and  time  constant  requirements.  The  tape  search  unit  generates  an 
output  TTL  comparable  signal  which  is  a logic  "1"  during  the  playback 
interval.  This  is  used  to  gate  the  analog-to-digital  conversion  of  the 
accelerometer  data. 

3.2.5  Time  Code  Recorder 


The  output  from  the  time  code  reproduction  amplifier  in  the  magnetic 
tape  playback  unit  is  routed  to  a time  code  reader.  The  time  code 
reader  is  designed  to  read  the  recorded  time  when  played  back  from  the 
magnetic  tape. 

The  standard  format  is  IRIG-B  and  includes  BCD  seconds,  minutes  and 
hours.  By  using  the  standard  TRIG  Format  B,  the  reader  is  usable  in 
conjunction  with  other  systems  or  projects.  The  unit  automatically  de- 
codes the  serial  input  code  and  provides  a front  panel  display  and 
digital  outputs.  A synchronization  mode,  selectable  from  the  front 
panel,  provides  noise  rejection.  This  feature  allows  the  selection  of 
"zero,"  "one,"  "two,"  "four,"  "eight"  or  infinite  time  frames  bypass  if 
input  errors  are  detected.  This  feature  is  very  important  when  digital 
outputs  are  used  or  tape  search  operation  is  desired. 
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3.2.6  Recorder 

The  magnetic  tape  recorder/reproducer  is  capable  of  handling  14 
channels  of  analog  data  at  seven  servo-controlled  tape  speeds  with 
magnetic  heads  that  are  TRIG  compatible.  The  methods  of  data  handling 
include  13  channels  of  FM  {H0%  deviation)  for  accelerometer  outputs, 
one  channel  of  direct  (analog)  for  recording  and  playing  back  time  code 
data  and  one-half  track  for  voice  annunciation.  The  main  recorder  case 
is  set  up  for  recording  only  FM  channels,  recording  and  playing  back 
voice  channels,  recording  only  time  codes  and  the  28-volt  dc  input  power 
supply  for  portable  aircraft  power  operation. 

In  an  operational  environment  this  recorder  unit  and  the  time  code 
generator  are  mounted  on  vibration  shock  mounts  for  easy  installation 
and  operation  on-board  the  aircraft.  The  magnetic  tape  recorder  is  run 
at  the  3-3/4  inch  per  second  speed  to  get  the  required  frequency  response 
out  to  1,000  Hz  for  the  data  analysis  phase.  This  speed  gives  the 
recording  system  the  capability  for  handling  about  3 hours  of  continuous 
data  per  tape. 

The  data  analysis  phase  of  the  system  utilizes  the  same  magnetic 
tape  recorder  unit  from  the  acquisition  phase  except  that  the  power 
supply  can  now  be  transformed  to  117-volt  ac  operation.  The  outputs  of 
each  channel  off  the  magnetic  tape  heads  are  routed  to  a portable  auxil- 
iary housing  that  has  its  own  117-volt  ac  power  supply.  This  housing 
has  all  the  reproduction  electronics  for  the  13  FM  channels  and  the  one- 
time code  channel  built  into  it. 

3.2.7  Time  Code  Generator 

The  time  code  generator  is  constructed  in  a standard  one-half  ATR 
configuration  per  Military  Specification  MS  91403  for  its  mounting  on 
tray  and  shock  mounts.  This  generator  is  capable  of  outputting  IRIG-B 
standard  time  (hours,  minutes  and  seconds)  for  correlating  accelerometer 
data.  This  unit  can  be  set  up  with  any  start  time  that  may  be  useful  to 
the  test  personnel  and  has  a visual  LED  display  on  its  front  panel  for 
ease  of  operation. 

3.2.8  Accelerometers  and  Amplifiers 

In  the  data  acquisition  phase  of  the  vibration  program,  measurements 
of  vibration  environments  will  be  accomplished  by  mounting  accelerometers 
in  aircraft  or  cargo  locations  of  concern  and  recording  the  acceleration 
time  histories  on  magnetic  tape.  The  system  initially  uses  two  low  mass 
triax  accelerometers  with  their  own  internal  preamplifiers  to  help 
eliminate  cable  noise  and  length  dependence.  These  accelerometers  have 
a battery  power  supply  and  external  amplifier  on  each  axis  to  provide  an 
external  X50  gain  before  inputting  to  the  recorder  unit.  This  external 
gain  is  needed  to  keep  low  acceleration  values  from  being  recorded 
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below  the  recorder's  noise  floor.  The  reason  for  this  is  that  the 
recorder  will  record  only  down  to  about  0.3  volt  in  the  FM  mode.  Any 
type  accelerometer  can  be  utilized  with  the  system. 

3.3  Software 

The  broad  programming  objectives  were  to  reformat,  by  the  use  of  a 
programmable  calculator,  the  outputs  from  an  analog  to  digital  converter 
to  a format  which  allows  the  project  engineer  a choice  in  the  area  of 
data  parmeter  presentation. 

' The  one  digital  technique  which  allows  the  system  to  operate  is  the 

use  of  an  algorithm  for  the  fast  computation  of  Fourier  transformations 
I on  a programmable  calculator.  The  algorithm,  usually  referred  to  as  a 

3 Fast  Fourier  Transform  (FFT)  procedure,  has  a great  impact  on  vibra- 

j:  tional  data  analysis  because  of  its  application  to  the  calculation  of 

I spectral  density  functions. 

Sampling  for  digital  data  analysis  is  usually  performed  at  equally 
spaced  intervals.  The  problem  then  is  to  determine  an  appropriate  sam- 
pling interval  h.  Consider  a continuous  record  which  is  sampled  such 
that  the  time  interval  between  sample  values  is  h seconds.  The  sampling 
rate  is  then  1/h  samples  per  second.  However,  a minimum  of  two  samples 
per  cycle  are  required  to  define  a frequency  component  in  the  original 
data.  Hence,  the  highest  frequency  which  can  be  defined  by  sampling  at 
a rate  of  1/h  samples  per  second  is  l/2h  Hz.  Frequencies  in  the  orig- 
inal data  above  l/2h  Hz  will  be  folded  back  into  the  frequency  range 
from  0 to  l/2h  Hz,  and  will  be  confused  with  data  in  the  lower  range. 

So  the  highest  frequency  that  can  be  used  for  a fixed  sampling  interval 
h is:  fc=l/2h  which  is  called  the  Nyquist  frequency  or  folding  frequency. 
In  this  system  there  is  a requirement  for  frequencies  up  to  1,000  Hz. 

So  h has  to  be  smaller  than  .0005  second.  On  the  other  hand,  sampling 
at  points  much  closer  will  yield  highly  redundant  data  and  thus  increase 
computer  time.  So  h was  chosen  to  be  .0004  second  for  easy  generation 
from  a 100  kHz  source.  The  tape  recorder,  because  of  the  3-1/2  inch  per 
second  speed,  will  cut  off  all  frequency  above  1.25  kHz  and  act  as  an 
antialiasing  filter  for  frequencies  above  this  Nyquist  frequency.  This 
means  that  the  analog-to-digital  converter  is  designed  to  sample  the  in- 
put analog  voltage  at  2.5  kHz. 

The  bipolar  two's  complement  code  from  the  12-bit  analog-to-digital 
! converter  is  read  into  the  calculator  using  a 16-bit  parallel  interface 

in  the  buffered  I/O  mode.  A fast  read/write  buffer  is  used  to  simply 
collect  the  data  as  fast  as  possible  without  spending  the  time  to  convert 
data  to  the  calculator's  internal  format.  This  work  is  done  later  by 
! the  input  program  after  the  data  has  been  collected.  This  buffered  I/O 

i mode  uses  an  8,000-word  buffer  to  read  about  three  seconds  of  digitized 


The  first  job  the  software  has  is  to  format  the  input  from  the 
analog-to-digital  12-bit  converter  to  a form  that  can  be  used  by  the 
Fast  Fourier  Transform  (FFT)  procedure  program.  The  converter  can 
handle  a voltage  swing  of  -10  volts  to  +10  volts  on  its  anlog  input  so 
its  digital  output  encodes  this  bipolar  voltage  into  a two's  complement 
code.  The  two's  complement  code  is  mathematically  achieved  for  negative 
numbers  or  voltages  by  complementing  corresponding  positive  number  and 
adding  a "1"  to  the  least  significant  bit.  The  resulting  set  is  indeed 
an  arithmetically  consistent  system  and  yet  has  only  one  zero  state. 

This  is  illustrated  in  Figure  6. 

When  using  the  16-bit  interface  in  the  fast  read/write  mode,  it 
operates  with  negati ve-true  logic  so  the  input  program  must  compensate 
for  inverted  logic. 

The  input  program  first  checks  the  status  line  from  the  A to  D to 
see  if  data  are  ready  to  be  written  into  the  8K  buffer.  After  this 
buffer  is  filled,  the  program  takes  512  words  from  it  and  places  these 
words  into  an  array.  Next  the  words  in  the  array  are  inverted  to  compen- 
sate for  inverted  logic  and  the  values  checked  to  see  if  greater  than 
2047.  If  greater,  the  number  is  negative  and  4096  is  subtracted  from 
it.  This  changes  the  number  range  to  -2048  to  +2047. 

Next,  the  numbers  are  divided  by  204,7,  which  changes  the  number 
range  again  to  the  new  range  of  -10  tc  +10.  This  range  corresponds  to 
the  voltages  seen  on  the  analog  input  to  the  A/D  converter.  Before  this 
has  happened,  the  calculator  has  asked  which  accelerometer  was  used  on 
this  channel  and  places  the  sensitivity  factor  (volts  per  g acceleration) 
into  variable  G.  This  value  is  used  to  convert  voltages  from  the  accel- 
erometers to  acceleration  to  be  used  by  the  Fast  Fourier  Transform 
program. 

The  input  program  stores  this  corrected  512-word  array  onto  tape 
and  takes  another  512  words  from  the  buffer  for  conversion.  After  the 
program  has  run  through  the  entire  buffer,  there  are  sixteen  512-word 
files  on  tape  for  the  FFT  program  to  operate  on. 

3.3.1  FFT  Procedural  Steps 

The  following  steps  were  used  tc  compute  power  spectra  estimates 
via  FFT  procedures.  Assume  that  the  sample  size  for  the  data  sequence 
Xn  is  initially  of  arbitrary  size  N which  in  our  case  is  512. 

a.  Truncate  the  data  sequence  or  add  zeros  so  that  N=2P  since 
Fast  Fourier  Transforms  operate  more  efficiently  in  this  matter  if  N is 
a power  of  two. 
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Ten  volts 

0111 

nil 

nil  = 2047 

full  positive 

0111 

nil 

• 

1110  » 2046 

-1  LSB 

Five  volts 

0100 

0000  ‘ 

0000  ■ 1024 

+half  full  scale 

0000 

0000  ’ 

0001  » 1 

+1  LSB 

Zero  volts 

0000 

0000 

0000  = 0 

Zero 

nil 

nil 

nil  = 4095 

-1  LSB 

Negative  five  volts 

1100 

0000  ’ 

0000  “ 3072 

-half  full  scale 

1000 

0000  ’ 

0001  “ 2049 

-FS+  LSB 

Negative  ten  volts 

1000 

0000 

0000  = 2048 

-FS 

Figure  6.  Two's  complement  12-bit  code 
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b.  Taper  the  resulting  sequence  using  the  cosine  taper  data  window 
of  Figure  7 to  take  into  account  the  effects  of  leakage.  Leakage  means 
that  components  of  f (t)  not  periodic  over  the  interval  NT  appear  to  leak 
out  into  portions  of  the  spectrum  adjacent  to  the  correct  spectral 
value.  This  is  caused  by  the  sharp  corners  of  the  rectangular  window 

w (t)  of  duration  NT  caused  by  the  digitizing  of  a discrete  number  of 
samples.  The  method  for  reducing  this  effect  would  be  simply  to  round 
off  the  sharp  corners  of  w (t)  by  using  the  cosine  taper  and  correcting 
for  this  window  after  doing  the  FFT. 

. n*N-l 

c.  Compute  Xl  where  Xu  = ^^^k» ' J ^ ^ Xp  exp  [-j  2Trkn] 

h n“o  N 

for  k=0,l,  N-1  where  X|(  are  the  Fourier  components  calculated  from  the 
FFT  procedure  using  the  adjusted  data  sequence  from  the  last  step. 

d.  Compute  the  power  spectrum  estimate  where  G|(  = (fk)  = 

where  the  frequency  of  each  estimate  is  f^  = ^ ^ and  K=0,1, 

...N-1,  N-1  is  the  number  of  the  estimate.  Note  that  the  arrays  are 
numbered  1 to  N,  not  0 to  N-1  in  the  actual  software. 

e.  Adjust  these  estimaj;es  by  a scale  factor  because  of  the  tapering 
used  in  step  2 by  replacing  Gk  by  (1/0.875)  Gk  if  the  cosine  tapering 
window  of  Figure  7 is  used. 

This  procedure  is  workable  for  only  a data  sequence  of  length  N. 

The  N is  the  number  of  points  used  by  the  FFT  algorithm  and  is  limited 
in  this  system  to  512  by  the  memory  size  of  the  calculator  for  storing 
arrays . 

The  method  used  to  get  around  this  is  to  divide  the  time  domain,  t, 
into  a number  of  segments  of  equal  length  N,  compute  a discrete  spectrum 
for  each  segment  using  the  FFT  algorithm  and  then  use  the  scaled  average 
of  the  squared  values  of  these  estimates  for  the  power  spectral  estimate. 
In  computing  the  discrete  spectrum  of  each  segment  a cosine  tapering  data 
window  was  impressed  over  each  segment  in  order  to  smooth  the  data 
sequence  for  reducing  leakage  as  described  before.  The  squared  Xk 
Fourier  components  for  each  segment  are  averaged  in  program  and  the 
final  spectral  estimate  put  on  tape  for  use  by  the  plotting  programs. 

In  this  method,  the  number  M of  frequency  bands  becomes  the  number 
of  independent  FFT  values  or  one-half  the  number  of  samples  in  a segment. 
The  number  of  segments  K was  chosen  to  be  16  because  this  is  the  maximum 
number  of  segments  the  calculator's  memory  and  tape  system  can  easily 
handle.  The  expression  for  the  power  spectral  estimate  now  becomes: 
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Assume  that  C^;  Is  defined  by  contlnous  function: 

T 


(t)  - COS 

2 


2 5jrt 


- cos^ 

T 


- 0 


T ^ < 2T 

5 - t - - 5 


2T  < 2T 
■5  ■ ' ■ 5 

21  ^ T 
5 " 2 


Otherwise 


Redefine  as  decrete  function: 


Where  N1  - [ (P2-P1)/10]+P1  and  N2-P2-[ (P2-P1)/10] 


C(N)  - sln2  52^) 


- 1 


. 2 .5?r(M-Pi) 
^Pi-Pi) 


PI  - N - m 


N1  - N - N2 


N2  - P2 


Otherwise 

FIGURE  7.  Cosine  taper  data  window 
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N-1 


*^k  = K = 0.1 

where  T = sampling  interval  in  seconds 

M = no.  frequency  values  from  0 to  1/2  Thz 

K = no.  segments 

W = scale  factor  for  adjusting  tapering 

Also  note  that  it  was  best  to  use  disjoint,  nonoverlapping  segments 
because  the  overall  record  length  of  f (t)  is  8,000  points  which  is  the 
maximum  size  buffer  the  calculator's  memory  can  handle.  The  output  from 
these  16  averaged  segments  are  recorded  on  a tape  file  for  use  by  plotting 
routines. 

Four  different  spectral  plots  are  used  in  this  system  to  analyze 
the  data;  these  are:  power  spectral  density,  acceleration,  double 
amplitude  and  double  amplitude  1og>log. 

a.  The  power  spectral  density  (PSD)  versus  frequency  program 
performs  no  additional  processing  on  the  data  contained  in  the  FFT 
output  file,  it  only  automates  the  job  of  plotting. 

b.  The  acceleration  versus  frequency  program  converts  PSD  data  to 
units  of  g rms  and  plots  that  data. 

c.  The  double  amplitude  versus  frequency  plot  program  integrates 
acceleration  values  derived  from  the  PSD  function  to  obtain  the  corre- 
sponding displacement  (in  inches)  at  each  frequency. 

d.  The  double  amplitude  log-log  uses  the  same  technique  as  the 
double  amplitude  versus  frequency  plot  program  except  in  log  scale  on 
both  axis  and  with  a MIL-STD-810C  vibration  curve  superimposed  on  it. 

The  curve  selected  is  curve  AW  for  equipment  category  6.  This  is  the 
curve  which  is  used  for  test  items  that  are  air  transported  by  the  Air 
Delivery  Engineering  Branch.  This  gives  the  project  engineer  a quick 
way  of  checking  if  all  the  data  points  in  this  segment  of  time  data 
falls  below  this  curve  and  passes  the  military  standard  it  is  being 
tested  against. 

4.0  RESULTS 

The  Airborne  Recording  System  has  successfully  utilized  the  versa- 
tility and  data  processing  capability  of  a modern  programmable  calculator 
for  vibration  analysis.  This  approach  makes  possible  the  addition  of 
other  analog  or  digital  tranducers  for  other  requirements  without  the 
procurement  of  different  analysis  instrumentation. 
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The  use  of  the  Fast  Fourier  Transfer  program  in  a calculator  this 
memory  size  has  limited  analysis  to  3.2  seconds  of  data  at  one  time. 

This  amounts  to  about  8,000  data  points  in  6 minutes  of  calculator's 
execution  time.  So  the  first  step  that  should  be  done  is  an  analog 
strip  recording  of  time  history  from  the  recorder  to  find  the  times  of 
interesting  events.  Then  these  times  should  be  used  for  detailed  analy- 
sis on  the  calculator. 

The  Airborne  Recording  System  was  first  used  during  the  air  transporta- 
tability  testing  of  the  Patriot  missile  which  utilized  a CH-47  Chinook 
helicopter.  The  acquired  spectral  vibrational  data  is  briefly  summarized 
in  the  plots  illustrated  in  Figures  8 through  11. 

The  Power  Spectral  Density  (PSD)  function  is  an  indication  of  the 
amount  of  vibrational  energy  at  a particular  frequency  at  a particular 
location  of  the  test  item  and  is  depicted  in  Figure  8.  The  test  item 
has  been  subjected  to  random  vibrations  along  each  of  three  mutually 
exclusive  axis  and  measured  by  triaxial  accelerometers,  placed  at  prede- 
termined locations.  The  Power  Spectral  Density  function  of  any  test 
item,  under  any  flight  environment,  shall  not  deviate  predetermined 
specifications  by  more  than  a predetermined  percentage  as  specified  in 
MIL-STD-810C.  In  addition,  the  PSD  function  provides  the  various  reso- 
nance peaks  that  the  test  item  is  being  subjected  to  aboard  the  aircraft. 

The  PSD  function  is  the  most  important  parameter  in  vibrational  analysis 
of  air  transportability  testing,  since  most  other  important  and  required 
parameters  can  be  derived  from  the  PSD  function. 

The  Acceleration  G RMS  function,  plotted  in  Figure  9,  is  used  to 
determine  resonance  response  peaks  during  a test  flight  in  air  transpor-  i 

tability  testing.  The  frequencies  at  which  significant  resonance  re-  I 

sponse  peaks  occur  are  identified  for  each  triaxial  accelerometer  mounted  j 

on  the  test  item.  This  information  is  then  utilized  for  resonance  dwell 
tests  to  determine  if  a test  item  can  withstand  maximum  dynamic  stress 
(endurance  tests)  as  defined  in  MIL-STD-810C.  | 

The  Double  Amplitude  function,  illustrated  in  Figure  10  is  required 
to  determine  the  maximum  movement  of  the  test  item  at  low  frequencies 
during  any  air  transportability  test.  This  measurement  is  required  for 
tie-down/stress  design  considerations. 

) 

The  Double  Amplitude  LOG  function,  plotted  in  Figure  11,  is  superim-  j 

posed  with  and  compared  to  a MIL-STD-810C  vibration  analysis  curve  for  ] 

the  air  transportability  equipment  category.  This  comparison  determines 
the  number  of  data  points  and  their  frequency  which  fall  within  the  MIL- 
STD-810C  curve  and/or  to  determine,  for  those  data  points  outside  the 
curve,  how  close  they  are  to  the  worst  case  curve. 
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Figure 


5.0  CONCLUSIONS 


The  Airborne  Recording  System  (ARS)  development  project  has  been 
successfully  completed  and  the  following  conclusions  arise: 

a.  The  development  of  the  ARS  instrumentation  system  provides  the 
project  engineer  with  a new  flexibility  in  the  area  of  spectral  data 
computation  and  presentation. 

b.  The  hardware  and  software  techniques  described  herein  provide 
an  alternate  method  for  augmenting  the  capabilities  of  existing  narrow- 
band,  swept  frequency  and  spectral  analysis. 

c.  All  of  the  instrumentation  procured  and  in-house  work  accom- 
plished satisfies  all  performance  requirements  made  by  the  Air  Delivery 
Engineering  Branch  of  YPG  in  terms  of  system  portability,  hardware 
reliability,  availability  and  maintainability  and  data  turn-around 
times. 


d.  Further  reduction  of  execution  turn-around  time  could  be  accom- 
plished with  a machine  language  (binary)  form  of  the  Fast  Fourier 
algorithm. 

6.0  RECOMMEDATIONS 

The  Airborne  Recording  System  (ARS)  has  been  successfully  developed 
within  the  technical  and  funding  constraints  as  specified  in  the  submitted 
proposals.  However,  a forecast  of  near  future  development  in  state-of- 
the-art  hardware  allows  us  to  make  the  following  recommendations: 

a.  Implement  a separate  developmental  task  to  expand  the  present 
Airborne  Recording  System  capabilities.  This  can  be  accomplished  through 
several  concurrent  subtasks.  The  execution  time  of  the  ARS  can  be 
speeded  up  with  a machine  language  (binary)  form  of  the  Fast  Fourier 
algorithm.  The  use  of  a floppy  disk  drive  for  recording  the  incoming 
data  in  larger  files  with  faster  access  time  will  reduce  the  time  and 
number  of  operator  actions  required  over  the  present  cassette  system. 
Adding  a quick-look  capability  by  utilizing  a video  monitor  will  cut 
down  time  required  for  printing  plots.  In  addition,  a high-speed  peak 
sample-hold  amplifier  w.’ll  cut  down  on  the  time  required  for  searching 
the  data  on  the  tape.  Tuis  hardware  will  provide  greater  system  flexi- 
bility, modularity  and  efficiency  to  facilitate  operator  interface  with 
input/output  devices. 
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b.  Modify  the  present  ARS  to  provide  a miniaturized  version  of  the 
airborne  components  of  the  ARS.  This  task  will  provide  a one-man  porta- 
ble ARS.  This  can  be  accomplished  with  a recorder/ reproduce  subsystem 
vhat  will  become  available  off-the-shelf  in  the  near  future.  The 
record/reproduce  subsystem  size  and  weight  will  be  reduced  by  a factor 
of  4.  Furthermore,  the  recent  development  of  FFT  modules  and  I/O  controller 
circuitry  will  allow  software  development  in  machine  language  for  perma- 
nent storage  on  read-only-memories  for  the  Data  Analysis  subsystem  of 
the  ARS. 


SECTION  II.  APPENDICES 
APPENDIX  A - ARS  OPERATING  PROCEDURES 

1.0  RECORDER/REPRODUCER  SYSTEM 

1.1  General  Preoperation  Procedures 

a.  Set  the  CALIBRATE  switch  (left  interface  panel)  to  OFF. 

b.  Set  the  TAPE/TONE  switch  (right  interface  panel)  to  TONE  for 
tone  wheel  signal  reference  of  the  capstan  servo  which  should  serve  most 
purposes. 


c.  Make  sure  the  tape  SPEED  switch  is  in  the  3-3/4  in. /sec  position. 

d.  Make  sure  SPEED  SELECTION  JACKS  for  equalization  SPEED  on  left 
side  has  a jumper  pin  in  the  jack  corresponding  to  3-3/4. 

e.  Make  sure  SPEED  SELECTION  JACKS  for  equalization  SPEED  B on 
right  side  has  a jumper  pin  in  the  jack  corresponding  to  15. 

f.  Mount  the  tape  reels  by  raising  the  locking  latch,  setting  the 
reel  over  the  hub  and  snapping  the  locking  latch  down. 

g.  Thread  the  tape  past  the  head  assemblies  and  around  the  capstan 
rollers  as  illustrated  on  the  bottom  of  the  head  cover  when  it  is  flipped 
up. 


1.2  Record  Operation  Procedures 

a.  After  checking  for  installation  of  28-volt  dc  power  supply, 
insert  record  transport  unit  into  its  mounting  tray.  See  recorder 
manual  for  instructions  in  the  change  of  the  power  supply  procedures. 

b.  Connect  coax  cable  from  each  accelerometer  to  an  input  BNC  con- 
nector on  interface  panels  and  label. 

c.  Connect  coax  cable  from  time  code  generator  to  channel  8 input 
BNC  connector  on  interface  panel. 

d.  Connect  power  cable  to  recorder  power  connector  and  other  end 
to  28-volt  dc  female  connector  next  to  recorder  tray. 

e.  Connect  headphone  and  mike  to  VA  input,  VA  output  and  voice  re- 
corder act  coax  BNC  connectors. 

f.  Connect  power  cable  from  the  time  code  generator  to  the  other 
28-volt  dc  female  connector  next  to  recorder  tray. 
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g.  Connect  power  cable  from  the  two  female  connectors  next  to  re- 
corder tra>  to  aircraft  power,  observing  proper  polarization  of  (-)  and 
(+). 

h.  Set  POWER  switch  on  time  code  generator  to  the  ON  position. 

i.  On  the  time  code  generator  the  nine-digit  display  shows  the 
accumulated  time  count.  The  DECADE  switch  selects  the  counter  decade 
(units-of-seconds,  units-of-minutes,  etc)  to  be  preset;  position  1 sel- 
ects the  units-of-seconds  decade  and  so  on  from  the  least  significant 
digit.  The  TIME  switch  selects  the  digit  to  be  entered  into  the  selected 
decade.  The~3ET  push  button  is  pressed  to  enter  the  selected  digit. 

The  START  push  button  is  used  to  manually  start  the  clock,  the  STOP  push 
button  is  used  to  manually  stop  the  clock.  The  RESET  push  button  is 
used  to  manually  reset  the  time  counters. 

j.  Turn  accelerometers  amplifier  boxes  ON  and  gain  switches  to 

X50. 

k.  To  start  the  recording  of  data,  set  POWER  switch  on  recorder  to 
ON.  Press  the  DRIVE  FORWARD  and  the  RECORD  push  buttons  at  the  same 
time. 


NOTE:  If  tape  does  not  begin  to  move,  check  that  tension  of  the 
tape  has  moved  the  tension  arms  away  from  their  extreme  released  position. 

1.3  Playback  Operation  Procedures 

a.  Check  for  installation  of  110-volt  ac  power  supply.  See  recorder 
manual  for  instructions  on  how  to  change  the  power  supplies. 

b.  Connect  the  coaxial  cables  as  required  for  channels  to  be 
reproduced;  the  cable  sequence  is  from  the  output  of  the  Model  5600 
recorder  to  the  input  of  the  auxiliary  housing  and  from  the  auxiliary 
housing  output  to  the  external  data  processing  equipment. 

c.  Connect  the  power  cables  from  the  auxiliary  housing  and  recorder 
to  facility  power. 

d.  Connect  the  auxiliary  cable  from  the  AUX  connector  on  the  Model 
5600  recorder  to  the  AUX  connector  on  the  auxiliary  housing. 

e.  Set  POWER  switch  on  both  units  to  ON. 

f.  See  tape  search  operation  procedures  for  rest  of  playback 
procedures. 
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2.0  TAPE  SEARCH  SYSTEM 


2.1  Turn  On  Procedure 

a.  Connect  the  interconnect  cable  from  J3  on  the  Model  9240  tape 
search  and  control  unit  to  the  Model  9200A  time  code  generator/ translator 
J4  connector. 

b.  Connect  tape  search  units  to  the  02  connector  tape  transport 
remote  operation  connector. 

c.  Place  the  standard  1 kHz  modulated  IRIG  B code  signal  from 
recorder/reproducer  time  code  channel  output  BNC  connector  to  BNC 
connector  01  on  Model  9200A  time  code  generator/translator. 

d.  Apply  115-volt  ac  power  to  06  connector  on  Model  9200A  time 
code  generator/translator. 

e.  Rotate  the  Mode  switch  to  the  TRANS  position  on  time  code 
genera tor/ translator. 

f.  Place  the  FWD/REV  switch  to  the  FWD  psition  on  time  code 
generator/translator. 

g.  Place  the  INPUT+/-  toggle  in  the  normal  (+)  position  on  the 
time  code  generator/transi ator . 

h.  Insure  that  the  FILTER  switch  is  in  the  position  on  the  time 
code  generator/translator. 

i.  Apply  115-volt  ac  power  to  J6  connector  on  the  tape  search  and 
control  unit. 

j.  The  Model  9240  tape  search  and  control  unit  and  the  Model  9200A 
time  code  generator/translator  are  powered  when  the  PWR  switch  is  placed 
to  the  ON  position  on  both  units. 

k.  Allow  approximately  one  minute  after  turn-on  for  AGC  to  settle; 
the  time  code  generator/translator  will  automatically  synchronize  to  the 
input  time  code. 

2.2  Manual  Search  Control  on  Tape  Search  System 

a.  Set  operational  mode  switch  to  MANUAL. 

b.  Between  each  manual  command,  a manual  stop  command  must  be 
initiated  by  pressing  STOP  on  the  tape  search  and  control  unit. 
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c.  STOP:  Press  STOP  to  stop  transport. 

d.  FWD:  Press  FWD  to  place  transport  in  FWD  mode. 

e.  FAST:  Press  FAST  pressing  either  FWD  or  REV  to  place  transport 
in  the  FAST  mode. 

2.3  Search  to  Start 

The  tape  search  unit  provides  for  time  increment  selection  from  the 
front  panel  by  means  of  thumbwheel  switches.  Thumbwheel  switches  (for 
hours,  minutes,  seconds)  are  provided  for  START  and  STOP  times.  Five 
operational  modes  are  selectable  from  front  panel. 

a.  Set  operational  mode  switch  to  the  SEARCH/START  position. 

b.  Preset  the  START  time  thumbwheel  switches  to  the  desired  times. 

c.  Press  START,  the  transport  will  stop  just  below  the  present 
start  time.  Once  stopped,  the  system  now  awaits  further  instructions. 

2.4  Single  Cycle  Mode 

The  tape  is  moved  ahead  of  start  time  and  is  stopped.  "Playback" 
is  automatically  activated  and  transport  is  placed  into  the  forward 
mode.  The  transport  is  stopped  at  the  detection  of  STOP  time  and  is 
ready  for  a new  operation.  This  is  the  normal  mode  for  reading  data  off 
of  tape  into  the  calculator. 

a.  Set  operational  mode  switch  to  the  SINGLE  CYCLE  position. 

b.  Preset  the  desired  START  time  on  the  START  TIME  thumbwheel 
switches  located  on  the  front  of  the  tape  search  and  control  unit. 

c.  Preset  the  desired  STOP  time  on  the  STOP  TIME  thumbwheel 
switches  located  on  the  front  of  the  tape  search  and  control  unit. 

d.  NOTE:  Four  seconds  difference  is  required  between  START  and 
STOP  time  to  fill  calculator  memory. 

2.5  Recycle  Mode 

This  mode  provides  for  repeated  single  cycle  operation.  The  proce- 
dures is  the  same  as  single  cycle  except  that  the  operational  mode 
switch  is  in  the  RECYCLE  position.  The  "STOP"  pushbutton  must  be  acti- 
vated to  stop  "RECYCLE"  operation. 
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3.0  CALCULATOR  SYSTEM 


3.1  Turn  ON  Procedure 

a.  Connect  the  power  line  Into  ac  power  outlets  for  printer  and 
calculator. 

b.  Switch  the  calculator  ON  using  the  switch  on  the  right-hand 
side  of  the  calculator. 

c.  Insert  the  tape  cartridge  so  that  the  label  on  the  cartridge 
faces  the  back  of  the  calculator. 

d.  Push  system  command  key  REWIND  to  rewind  the  tape  cartridge  to 
its  beginning. 

e.  Connect  the  HIB-bus  interface  cable  to  the  I/O  connector  at  the 
rear  of  the  printer,  and  to  the  calculator  rear  peripheral  connector. 

f.  Switch  the  OFF/ON  switch  located  on  the  right  front  of  the 
printer  to  the  ON  position. 

g.  Connect  the  16-bit  98032A  interface  cable  from  the  digitizer  to 
one  of  the  calculator  rear  peripheral  connectors. 

h.  Connect  the  power  cord  from  the  digitizer  into  an  ac  power 
outlet. 

3.2  Calculator  Operation 

a.  Type  idf  0 and  press  execute  key;  this  loads  the  input  program 
into  the  calculator's  memory. 

b.  Connect  coax  cable  from  the  desired  channel  of  the  recorder/ 
reproducer  to  the  digitizer's  BNC  connector. 

c.  Set  the  tape  search  system  for  desired  START  and  STOP  times. 
Note  that  it  must  be  four  seconds  difference. 

d.  Push  RUN  on  calculator  keyboard. 

e.  Start  tape  search  system  operation  as  described  before. 

f.  The  calculator's  display  will  ask  what  accelerometer  number  was 
utilized  on  this  channel.  Type  in  the  number  from  the  keyboard  and  push 
the  CONTINUE  key. 
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g.  The  calculator's  display  will  ask  the  gain  used  with  that 
accelerometer.  "1"  means  times  1 and  "2"  means  times  10.  Type  in  the 
number  from  the  keyboard  and  push  the  CONTINUE  key. 

h.  When  the  "lazy  T"  shows  on  the  display  the  input  program  is 

done. 

i.  Push  the  ERASE,  the  A and  the  EXECUTE  keys  sequentially  on  the 
calculator's  keyboard.  This  erases  the  entire  calculator  memory. 

j.  Type  idf  1 and  press  EXECUTE  key;  this  loads  the  FFT  program 
into  the  calculator's  memory. 

k.  Push  RUN  key  on  the  calculator's  keyboard  to  start  the  program. 

l.  When  the  "lazy  T"  shows  on  the  display,  the  FFT  program  is 

done. 

m.  Push  the  ERASE,  the  A and  the  EXECUTE  keys  sequentially  on  the 
calculator's  keyboard.  This  again  erases  the  entire  calculator  memory. 

n.  Type  Idf  N where  N=19  for  power  spectral  density  plot,  N=20  for 
accleration  plot,  N=21  for  double  amplitude  plot  and  N=22  for  double 

I amplitude  log  plot.  Press  EXECUTE  key  for  loading  desired  plot  program 

' into  the  calculator's  memory. 

i 

[ 0.  Push  RUN  key  on  the  calculator's  keyboard  to  start  the  plotting 

program. 

p.  When  plotting  is  done  the  calculator's  display  will  ask  if 
another  copy  of  the  plot  is  desired.  Type  "1"  for  YES  and  "0"  for  NO 
and  press  the  CONTINUE  key. 

q.  When  the  "lazy  T"  shows  on  the  display  the  plotting  program  is 

done. 

r.  Push  the  ERASE,  the  A and  the  EXECUTE  keys  sequentially  on  the 
calculator's  keyboard.  This  again  erases  the  entire  calculator  memory. 

s.  Repeat  steps  n.  through  r.  for  all  different  plots  desired  for 
that  data  segment. 

t.  Repeat  all  operation  steps  for  each  new  data  channel  or  segment 
being  analyzed. 
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u.  Check  for  paper  in  printer.  To  load  paper,  first  move  paper 
release  lever  back,  place  the  paper  behind  the  platen  with  the  left  edge 
of  paper  at  the  desired  position.  Lift  the  paper  bail  up  then  turn  the 
platen  knob  clockwise  to  feed  the  paper  around  the  platen.  Push  the 
paper  bail  down  and  ensure  that  the  paper  travels  behind  the  bail.  To 
remove  or  reposition  the  paper  in  the  printer,  move  the  paper  release 
lever  -forward.  With  the  lever  forward  the  paper  is  free  to  move.  When 
the  paper  is  properly  positioned  (horizontally  and  verically),  move  the 
lever  back  for  details,  consult  the  printer  manual. 
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APPENDIX  B.  PROGRAM  LISTINGS 


( 


I 

COPY  AVAIIME  T!!  ""  ‘ ^ 

f ■*!  *’  ■ * y' 

. *i|;|jrDT'"^[}/(f5l(‘pROGRAM 

0:  dim  X[5121 ,A(61 
1;  buf  "hold", 8192, 2 
2:  ioEl3>L;if  l>0;jmp  0 
3:  tfr  13, "hold", 8192 
4:  rdsC'hold"  ) ♦K;  if  K*-l;jmp  0 
5:  .00975-A[l] 

6:  .00994*A[2] 

7:  . 00988-A[3] 

3:  .01015-A[4] 

9;  .0098-A[5] 

10:  .01078*A[6] 

11:  dsp  "EdTER  ACC.  TRANSDUCER  NUMBER" ;wait  1000 
12:  ent  R 
13:  AlRl*3 

14:  dsD  "ENTER  1 FOR  IX  SENSITIVITY  OR  2 FOR  lOX" ;wait  1000 
15:  ent  Q 
15:  if  C*l;gto  +3 
17:  10-1 
18:  g t o + 2 
19:  1-1 
20:  !*3*G 

21:  for  J=1  to  16 
22:  for  1=1  to  512 
23:  rdb("hold")*X[I] 

24:  -X[I]>X[I] 

25:  gto  +l;if  XI I ] <2048 ;gto  +2 
26:  X[11-4096*X(I] 

27:  next  1 
23:  for  1=1  to  512 
29:  X[I)/(G*204.7)*X[I] 

30:  next  I 
31:  1-A 
32:  512^B 
33:  {3-A)/10+A*C 
34:  3-(a-A)/10*E 
35:  for  I=A  to  B 
36:  gto  +l;if  I<»C;gto  +4 
37:  gto  +l;if  I>=E;gto  +3 
38:  wait  1 
39:  gto  +3 

40:  rad;sin(5*TT*(I-A)/{B-A))^F 
41:  F‘2*X[I]>XlI] 

42:  next  I 
43:  ref  J+1,X[*] 

44:  next  J 
45:  end 
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FAST  FOURIER  TRANSFORM  PROGRAM 


0:  512+N 

1:  dim  X [ 512 ] , Y [ 512] ,C [ 129] , D [2 ,256] 
2:  for  A=1  to  16 
3:  Idf  A+1,X[*] 

4:  N/4*G;3+1*H 

5;  H+1*J;G+J*K 

6:  int(log(N)/log(2)+.5)*L 

7:  2*it/N-»M 

8:  for  0=1  to  H 

9:  cos(  (0-1)  *M)  *0  [0] 

10:  next  0 

11;  for  0=1  to  L 

12;  2''(L-0)-P 

13;  2*P-*Q;N/Q+R 

14;  for  S=1  to  P 

15:  (S-1)*R+1-T 

16;  if  T<=H;gto  +4 

17;  -G[K-Tl+M 

18;  -C[T-G]+X 

19;  gto  +3 

20;  ClT]*M 

21;  -C[J-T]*X 

22;  for  U=0  to  N-Q  by  Q 

23;  U+S*V 

24:  V+PfY 

25;  XIV]-X[Y]+Z 

26:  Y(V]-Y[Y]*D 

27;  X[V]+X [Y]>XfV] 

28:  Y[V] +Y  [Y]-*-Y[V] 

29;  I1*Z-X*D-»XIY1 
30;  i4*D+X*Z+Y  [Y] 

31;  next  U 

32:  next  S 

33:  next  0 

34;  L+1-0 

35;  for  P=1  to  8 

36;  1-C(P] 

37;  if  P>L;gto  +2 
38:  2''(0-P)-CtP] 

39:  next  P 
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BESI  AVAIIME  COPY 

j 

i 

40 

1-0 

41 

for  0=1  to  C[3] 

42 

for  ^=0  to  C(7]  by  C(3] 

43 

for  C=P  to  C[5]  by  C(7] 

44 

for  R=Q  to  C[5)  by  C(6] 

45 

for  S=R  to  C[41  by  C151 

45 

for  T=S  to  C(3]  by  C[4] 

47 

for  U-I  to  C[2]  by  C[31 

48 

for  V=U  to  C[l]  by  C[2] 

49 

if  D>V;gto  +7 

50 

X [ D]  VI 

51 

X[\/]  ♦X[D] 

52 

M-^XW] 

53 

y[D]  VI 

54 

Y[y ] ♦y [D] 

55 

•1-y  [y/] 

55 

D+1*D 

57 

next  V 

58 

next  U 

59 

next  T 

60 

next  S 

61 

next  R 

62 

next  Q 

63 

next  P 

64 

next  0 

65 

0*3 

66 

.0004-*W 

67 

for  1*1  to  256 

68 

X(I] “2+y (I ]“2-X[I] 

69 

W*X[I ]/(N*.875)+D [1,1] ♦D[ 1,1] 

70 

next  I 

71 

ina  y 

72 

ort  A 

73 

next  A 

74 

for  1*1  to  256 

75 

I/(N*W) [2,1 ] 

76 

D[1,I)/16*D[1,I] 

77 

next  I 

78 

ref  18, D[*) 

79 

end 
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2:  0-C;1000-E;n-E;ldf  1R,D[*] 

3:  for  1=1  to  205 

4;  if  0[l,I)<3;gto  +2 

5:  D[l,I]-8 

5:  next  I 

7:  100  + F ;1000'*D 

3;  cll  'forTi' (14,11,11) 

9:  cll  'psiz' (8,9,1,2) 

10;  cll  'scl' (C-.2*D,E,-. 1*8,1. 1*3) 

11:  cll  'xaxis' (0, F,C,E) 

12:  cll  'yaxis' (C, . 1*3,0 ,1. 13) 

13:  for  G=C  to  E by  F 
14:  cll  'move'(G,0) 

15:  cll  'soace' (-3) 

16:  cll  'skio'(l) 

17;  if  G>E;qto  +? 

18 : fmt  2 , f 4 . 0 
19:  wrt  701. 2, G 
20;  next  G 

21:  cll  'move' ( .4*0, -.05*3) 

22:  wrt  70 1 , " FPEQUENCY-HEETZ  " 

23:  for  G=0  to  1.2*3  by  .1*3 
24:  cll  'move' (C-. 15*0  ,G) 

25:  fmt  4 ,e8. 1 
26;  wrt  701. 4, G 
27:  next  G 

28:  cll  'move' ( .05*0 ,1.1*3) 

29:  wrt  701, "POWER  SPECTRAL  DEU-G*G/nZ" 
30;  cll  'move'(C,0) 

31:  for  J=1  to  256 
32:  if  0[2, J]>=C;gto  +2 
33:  gto  +3 

34;  cll  'fnlt'(0[2,J],D[l,Jl) 

35:  if  D[2, J]>=E;gto  +2 
36:  next  J 

37:  dEO  "00  YOU  WANT  ANOTHER  PLOT" 

38:  wait  500 
39:  ent  W 
40 ; wtb  rO  ,13 
41:  if  7=l;gto  10 
42:  end 
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43:  "^ove”: 

44:  wtt  rO  ,27,65  , int  ( (nl-X)a/<54)  , int  ( ( r] -x)  U ) , i nt  ( (n3-Y)  V/G  4) 

, int{ (r2-Y)V) 

45:  ret 
46:  "nit": 

47:  wtb  rO  ,27,65  , int  ( (nl-X)U/64)  , int  ( ( cl-X)  J)  , irt  ( (-)2-Y)  V/64) 

, ini  ( ( n2-Y)\/) 

4 ' : if  n3=0 ; 45*n3 

49:  if  o3  = 4'5;wtb  r 0 , 2 7 , 3 2 ,0 , 0 , 0 , 6 

50:  wtb  t0,n3;wtb  cO  ,3 

51.  if  n3  = 4G;wtb  r 0 , 2 7 , 3 2 , 0 , 0 , 6 3 ,- 6 

5L : ret 

53:  "felt": 

5'.:  wtb  rO  ,27,97  , int  ( (ol-X)U/G4)  , int  ( (nl-X)U)  , int  ( (n2-Y)  v/64) 

, int { (n2-Y)V) 

55:  if  D3=0;46*n3 

55:  if  n3  = 4 6;v7tb  rO  , 2 7 , 3 2 ,0 , 0 ,0 , 6 

57:  wtb  r0,n3;wtb  r0,3 

53:  if  n3=46;wtb  rO ,27 ,3 2 ,0 ,0  ,5 3 ,-6 

55:  ret 

50:  "nsiz": 

51:  nl+[l;  n2^v7 

52:  wtb  r0,27,79,  int  (o4  * .1  20/6  4)  , nl  * 1 20  , int  (n3  * <^6 /64  ) ,r.3*Q6 
63:  ret 
64:  "pel": 

65:  120J/(n2-pl) ♦U 
66  : 95il/(  p4-r)3)  ♦V 
57:  o1+X;p3^Y 
63:  ret 
69:  "xaxis": 

70:  wtb  rO, 27, 46, 95, 0,5,9 

71:  if  p3=0  and  n4  = 0 j X*p3  ; X+120  T/U*r>4 

72:  if  p2=0 ; n4-p3>p2 

73:  wtb  rO ,27,65,  int  ( (p3-X) U/54)  , int ( (n3-x) U) , int ( (Pl-Y) V/64) 
,int((pl-Y)V) 

74:  P3>p5;wtb  r0,43;wth  rO ,3 

75:  wtb  rO  ,27 , 114  , int  (p2U/64)  , int  { ’"2U)  ,0  ,0  : wtb  rO  , 4 3 , 8 ; i pp  (p 
5+r2-*-r5)  >=p4 
76:  ret 
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77;  "Yaxis": 

73:  v.'tb  rO, 27, 46, 124, 0,3,0 

79:  if  n3  = 0 and  t4  = 0 ; Y+p3  ; Y+96d/Y-*’n4 

80:  if  o2=0 ; n4-o3+n2 

31;  wtb  rO  ,27,65  , int ( (ol-X)U/54)  , int ( (pl-X) U) , int ( (d3-Y)V/64) 

, int ( ( p3-Y)V) 

32;  n3 -►05  ; v^?tb  rO,43jv/tb  r0,3 

33;  wtb  rO , 27 , 114 , 0 , 0 , int ( d2 V/54) , int { p2V ) ; wtb  rO , 4 3 , 8 ; jmc  (p 

5-KC2  •►p5 ) > = o4 

34;  ret 

35:  "snace": 

35:  if  nl<0;qtc  +2 

37:  v;tb  rO,32;jrRp  2 ( (Dl-l-»pl)  =0) 

33:  wtb  rO,8;jiT'o  (Dl  + l-»nl)=0 
89:  ret 
90:  "skio": 

91 : if  oKO  ; q tc  +2 

92:  wtb  r0,10;jq'o  2 ( ( ol-l -►ol)  =0) 

93:  wtb  rO  ,27  , io  ; imo  (ol-i-l-^ol)  =0 
94:  ret 
95;  "fcro": 

96;  wtb  ,21 ,11 

97:  wtb  r0,27,84 

98:  if  nl=n  ; 13. 2-^ol ; Il-»n2-^n3 

99;  wtb  rO  ,27,37,  int  (120*nl/54)  ,120*pl 

100;  wtb  r0,27,75,int(96*n2/64) ,96*c2 

101:  wtb  r0,27,70,int(96*p3/64) ,95*03 

102:  ret 
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ACCELERATION  PLOT  PROGRAM 


0:  dim  D[2,256] 

1;  0-B;ldf  18, Dl*] 

2:  4.88+R 
3:  for  1=1  to  256 
4:  /(R*D[1,I ] )-L 
5:  if  L<B;gto  +2 
6 : L-»B 
7:  next  I 
8 ; pr t B 
9:  701-r0 

10;  0+C; 1000*E;1000*D;100-F 
11;  oil  'form' (14,11,11) 

12;  oil  'psiz' (8,9,1,2) 

13;  oil  'scl' (C-.2*D,E, -.1*8,1. 05*B) 
14;  oil  'xaxis' ( 0 , F,C,E ) 

15;  cll  'yaxis' (C, . 1*B,0 ,1B) 

16;  for  G=C  to  E by  F 
17;  cll  'move'(G,0) 

18;  cll  'space' (-3) 

19;  cll  'skip'(l) 

20;  if  G>E;gto  +3 
21;  fmt  2,f4.0 
22;  wrt  701. 2, G 
23;  next  G 

24;  cll  'move' ( .4*D,-. 05*B) 

25;  wrt  701, "FREQUENCY-HERTZ" 

26;  for  G*0  to  B by  .1*B 
27;  cll  'move' (C-.16*D,G) 

28;  fmt  4,e8.1 
29;  wrt  701 . 4 ,G 
30;  next  G 

31;  cll  'move' ( .05*0 ,B) 

32;  wrt  701 , "ACCELERATION  G RMS" 

33;  cll  'move'(C,0) 

34;  for  J=1  to  256 
35;  if  D[2, J)>=C;gto  +2 
36;  gto  +3 

37;  cll  'fplt' (D [2, J] , V(R*D[1,J] ) ) 
38;  if  D[2, J]>=E;gto  +2 
39;  next  J 

40;  dsp  "DO  YOU  WANT  ANOTHER  PLOT" 

41;  wait  500 

42;  ent  W 

43;  wtb  rO  ,13 

44;  if  W=l;gto  11 

45;  end 
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465  "move" 5 

47;  wtb  r0,27,65,int{(ol-X)U/64) ,int((pl-X)U) ,int((p2-y)V/64) 
,int((p2-Y)V) 

48:  ret 
49!  "pit"; 

50;  wtb  r0,27,65,itit((ol-X)U/64)  ,int({pl-X)U),int((D2-Y)V/64) 
,int({p2-Y)V) 

51:  if  p3»0;46*p3 

52:  if  p3»46;wtb  rO, 27 ,8 2 ,0 ,0 ,0 ,6 

53:  wtb  r0,p3!Wtb  rO ,8 

54:  if  p3'*46;wtb  rO , 2 7 ,8  2 ,0 ,0  ,6  3 ,-6 

55:  ret 

56:  "fplt": 

57:  wtb  r0,27,97,int((ol-X)U/64) ,int{(pl-X)U) ,int( (o2-Y)V/64) 
,int((p2-Y)V) 


58; 

if  p3«0;46*p3 

59: 

if  p3»46!wtb  rO ,27,8 2 ,0 ,0  ,0 ,6 

60: 

wtb  r0,p3»wtb  r0,8 

61: 

if  p3»46;wtb  rO ,27 ,8 2 ,0 ,0  ,6 3 ,-6 

62; 

ret 

63: 

"psiz": 

64: 

pi  *lif  p2^W 

65: 

wtb  r0,27,79, int (p4* 120/64) , p4*120 , int (p3*96/64) ,p3*96 

66: 

ret 

67: 

"scl": 

68: 

120W/(p2-pl)*U 

69: 

96H/(p4-p3)*V 

70: 

pi  ♦X  j p3  ♦Y 

71: 

ret 

72: 

"xaxia" ; 

73: 

wtb  rO, 27, 46, 95, 0,5, 9 

74: 

if  p3«0  and  p4»0>X*p3 »X+120W/U^p4 

75: 

if  p2«0; p4-p3*p2 

76: 

wtb  r0,27,65,int((p3-X)U/64) , int ( ( p3-X) U) , 

, int( (ol-Y)V/64) 

,int((pl-Y)V) 

77:  p3*p5;wtb  r0,43;wtb  r0,8 

78:  wtb  rO ,27,114,int(p2U/64) ,int{p2U) ,0,0;wtb  r0,43,8;jmp  (p 

5+p2*p5)  >*p4 

79:  ret 

30:  "yaxis"; 

81:  wtb  rO, 27, 46, 124,0 ,3,0 

82:  if  p3-0  and  p4-0! Y*p3 ; Y+96H/V*p4 
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8 3:  if  p2"0jp4-t»3*p2 

84*  wtb  rO,27,6S*lnt«t»l-X)U/54)  ,int{(ol-X)U)  ,int<(o3-y)V/64) 
,lnt((p3-Y)V) 

85:  p3*o5;wtb  cO,43|vrtb  r0,8 

86:  wtb  r0,27,114,0,0^int(p2V/64) ,int(p2V)jwtb  r0,43,3;jmp  (p 

5+p2*p5) >»p4 

87*  ret 

88:  •epace*! 

39:  if  pKOigto  ^2 

90:  wtb  r0«32|jinp  2 ( (pl-l^pl)  *0) 

91*  wtb  cO^SijmD  (oi+i^pl)"0 
92*  ret 
93:  "skip"* 

94*  if  pl<0}gto  >2 

95*  wtb  r0,l0i3inp  2 ( (pl-l*pl)  »0) 

96*  wtb  r0,27*10j3fBp  <pl4-l^pl)  *0 
97:  ret 
98:  "form"* 

99*  wtb  rO, 27*77 

100*  wtb  r0*27^84 

101:  if  Pl«0|l3*2*pl|ll*p2^p3 

102*  wtb  r0,27,a7,int(l20*pl/64) ,120*pl 

103*  wtb  r0*27,76,lnt(96"p2/64)  ,96*p2 

104*  wtb  r0,27,70,int(96*p3/84) *96*p3 

105*  ret 


1 


best  available  copv 


DOUBLE  AMPLITUDE  PLOT  PROGRAM 

0:  dim  D[2,256] 

1:  701-r0 

2;  0-K:;1000*E;  1000*D;4.88*R 
3:  384.6*/(2*R)/(2*n"2)+M 
4:  0-Q;ldf  18, D [*] 

5:  for  1=1  to  256 

6:  M*/D  [1,I]/D{2,I]“2*L 

7;  if  L<=Q;gto  +2 

8:  L+Q 

9:  next  I 

10:  O-^N 

11:  if  Q>=l;gto  +5 

12:  N+1*N 

13:  Q*10*N*P 

14:  if  int(P)#0;gto  +9 

15:  gto  -4 

16:  if  Q>10;gto  +3 

17:  Q-P 

18:  gto  +5 

19:  N-l+N 

20:  Q*10"N*P 

21:  if  P<10;gto  +2 

22:  gto  -3 

2 3:  (int(P)  + (int(P)#P)  )*10“(-N)  -B 
24:  100*F 

25:  oil  'form'(14,ll,ll) 

26:  oil  'psiz'(8,9,l,2) 

27:  oil  'scl' (C-.2*D,E,-.1*B,1.05*B) 

28:  oil  'xaxis' (0,F,C,E) 

29:  oil  'yaxis' (C,.1*B,0,B) 

30:  for  G=C  to  E by  F 
31:  oil  'move'(G,0) 

32:  oil  'space' (-3) 

33:  cll  'skip'(l) 

34:  if  G>E;gto  +3 
35:  fmt  2,f4.0 
36:  wrt  701. 2, G 
37:  next  G 

38:  cll  'move' ( .4*0,-. 05*B) 

39:  wrt  701, "FREQUENCY-HERTZ" 

40:  for  G=0  to  B by  .1*B 
41:  cll  'move' (C-. 16*0, G) 

42:  fmt  4,e8.1 
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43:  wrt  701. 4, G 
44:  next  G 

45:  cll  'move' ( .05*D,B) 

46:  wrt  701, “DOUBLE  AMPLITUDE-INCHES" 

47:  cll  'move'(C,0) 

48:  for  J*1  to  256 
49:  if  D(2,^]>»C;gto  +2 
50:  gto  +3 

51:  cll  'fplt'(Dl2,Jl,M*/D(l,J]A)l2,Jl*2) 

52:  if  D[2, J]>»Efgto  +2 
53:  next  J 

54:  dsp  "DO  VOU  WANT  ANOTHER  PLOT" 

55:  wait  500 
56:  ent  W 
57:  wtb  r0,13 
53:  if  W=l;gto  25 
59:  end 
60:  " move " : 

61:  wtb  r0,27,65,int((ol-X)U/64) ,int((pl-X)U) ,int((p2-Y)V/64) 
,int((p2-Y)V) 

62:  ret 
63:  "pit": 

64:  wtb  r0,27,65,int((pl-X)U/64) ,int((pl-X)U) ,int((p2-Y)V/64) 
,int((p2-Y)V) 

65:  if  p3*0;46-»p3 

66;  if  p3»46;wtb  r0,27,8 2,0 ,0 ,0 ,6 

67:  wtb  r0,p3;wtb  r0,8 

68;  if  p3»46jwtb  r0,27 ,8 2,0 ,0 ,6 3 ,-6 

69;  ret 

70:  "fplt": 

71;  wtb  r0,27,97,int((pl-X)U/64) ,int{(pl-X)U) ,int((p2-Y)V/64) 
,int((p2-Y)V) 

72:  if  p3=0;46*p3 

73:  if  p3“46;wtb  r0,27,82,0 ,0 ,0,i6 

74:  wtb  r0,p3;wtb  rO ,8 

75:  if  p3=46;wtb  rO, 27,82,0,0,63,-^6 

76;  ret 

77;  "psiz"; 

78:  pi*H;p2*W 

79;  wtb  r0,27,79,int(p4*120/64) ,p4*120 ,int(p3*96/64) ,p3*96 

80:  ret 

81*  " qcl" t 

82:  120W/(p2-pl)*U 

83;  96H/(p4-p3)*V 
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84:  pl*X;p3-»y 

85:  ret 

86:  "xaxis": 

87:  wtb  rO, 27, 46, 95, 0,5, 9 

88:  if  p3=0  and  p4=0; X*p3 ;X+120W/0+p4 

89:  if  p2=0 ; p4-p3+p2 

90:  wtb  t0,27,65,int({p3-X)U/64) ,int((p3-X)U) ,int( (pl-Y)V/64) 
,int{(pl-Y)V) 

91:  p3-»p5;wtb  r0,43;wtb  r0,8 

92:  wtb  rO ,27 ,114 , int (p2U/64) , int ( p2U) ,0 ,0 ; wtb  r0,43,8;jmp  (p 

5+p2*p5) >=p4 

93:  ret 

94:  "yaxis": 

95:  wtb  rO, 27, 46, 124,0 ,3,0 

96:  if  p3=0  and  p4=0; Y^p3 ; Y+96H/V*p4 

97:  if  p2=0 ; p4-d3+p2 

98:  wtb  r0,27,65,int((pl-X)U/64) ,int((pl-X)U),int({D3-Y)V/64) 
,int((p3-Y)V) 

99:  p3-*’p5;wtb  r0,43;wtb  r0,8 

100:  wtb  rO ,27 ,114 ,0 ,0 , int (p2V/64)  , int (p2V) ;wtb  rO ,4 3 ,8 ; jmo  ( 
p5+p2+p5) >=p4 
101:  ret 
102:  "space": 

103:  if  pl<0;gto  +2 

104:  wtb  r0,32;jmp  2 ( (pl-l+pl)  =0) 

105:  wtb  r0,8;jmp  (pl+l+pl)=0 
106:  ret 
107:  "skip": 

108:  if  pl<0;gto  +2 

109:  wtb  r0,10;jmp  2( (pl-l^pl)  =0) 

110:  wtb  r0,27,10;jmp  (pl+l^pl) =0 
111:  ret 
112:  "form": 

113:  wtb  r0,27,77 

114:  wtb  r0,27,84 

115:  if  Pl=0;13. 2-pl;ll-p2-p3 

116:  wtb  r0,27,87 ,int(120*pl/64) ,120*pl 

117:  wtb  r0,27,76,int{96*p2/64) ,96*p2 

118:  wtb  r0,27,70,int{96*p3/64)  ,96*p3 

119:  ret 


BBI  AVWIABIE  COPY 

DOUBLE  AMPLITUDE  LOG-LOG  PLOT  PROGRAM 


0 : dim  0(2,256] 

1:  701-r0 
2:  1*C 

3:  1000-*-E;4.88*R 
4;  384.6*/(2*R) /(2*tt'2)-M 
5:  Idf  18, D[*] 

6:  cll  'form' (14,11,11) 

7;  cll  'psiz ' ( 8 ,9 , 1,2 ) 

8;  cll  'scl'(-l,5,-6,l) 

9:  cll  'xaxis' (-5,1,0,4) 

10:  cll  'yaxis' (0,1, -5  ,0) 

11:  for  G=0  to  4 by  1 
12:  cll  'move' (G, -5) 

13:  cll  'sDace'(-3) 

14:  cll  'skip'(l) 

15:  fmt  2,f4.0 
16:  wrt  701. 2, G 
17:  next  G 

18:  cll  'move' ( 1. 5,-5. 3) 

19:  wrt  701, "LOG  FREQUENCY-HERTZ" 

20:  for  G=-5  to  0 
21:  cll  'move' (-.2, G) 

22:  fmt  4,f2.0 
2 3:  wrt  701 . 4 ,G 
24:  next  G 
25;  cll  'move' (2,0) 

26:  wrt  701, "DOUBLE  AMPLITUDE  LOG-INCHES" 

27:  for  J=1  to  256 

28;  if  D[2, J] <C;Tto  +4 

29;  M*/D(l,J]/Dl2,Jr2-*-N 

30:  cll  'olt' (loq(D(2,J] ) ,log(N) ) 

31;  if  D[2,J)>F;citc  +2 
32;  next  J 

33;  cll  'move' ( log ( 5 . 5) , log ( le-5) ) 

34;  cll  'folt' (log(5. 5) , loq( leO)  ,46) 

35:  cll  'char ' (46,2,6) 

36;  cll  'folt' ( log  ( 38)  , log  ( 2e-2)  ,46) 

37:  cll  'fnlt'(log(200) ,log(7.33432e-4) ) 
38:  cll  'folt'(log(200) ,log(le-5)) 

39;  cll  'move' ( log(33) ,log(8e-2) ) 

40;  wrt  701,"1.5G" 

41:  cll  'move' ( lon(500) , log(le-2)  ) 

42:  wrt  701 , "MI L-STD-8 lOC-CURVE  AW" 

43:  cll  'move' ( loa{ 500) , log(le-3)  1 
44;  wrt  701, "FOr  EQUIPMENT  CATEGORY  0" 

45;  dsp  "DO  YOU  WANT  ANOTHER  PLOT" 
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46:  wait  500 
47:  ent 
43:  wtb  r0,13 
49:  if  W=l;gto  6 
50:  end 
51:  "move": 

52:  wtb  r0,27,65,int((pl-X)U/64) ,int((pl-X)U) ,int((D2-Y)V/64) 
, int( (p2-Y)V) 

53:  ret 
54:  "olt": 

55:  wtb  r0,27 ,65 , int { (d1-X)U/64)  , int ( (pl-X)U) , int ( (p2-Y)V/64) 
, int( (p2-Y)V) 

56:  if  p3=0;46*d3 

57:  if  p3=46;wtb  rO , 2 7 , 8 2 ,0 , 0 , 0 , 6 

58:  wtb  r0,p3;wtb  r0,8 

59:  if  p3=46;wtb  rO  ,2 7 , 8 2 ,0 , 0 ,6 3 ,-6 

60:  ret 

61:  "fplt": 

62:  wtb  rO ,27,97 , int ( (d1-X)U/64) , int( (pl-X)U) , int( (d2-Y)V/64) 
, int( (p2-y)V) 

63:  if  p3=0;46+p3 

64:  if  p3=46;wtb  rO , 2 7 , 8 2 ,0 , 0 , 0 , 6 

65:  wtb  r0,D3;wtb  r0,8 

66:  if  p3=46;wtb  r 0 ,27 ,8  2 ,0 , 0 ,6  3 ,-6 

67:  ret 

68:  "char": 

69:  if  p2  = 0 ; 5^p2 ; 0*d3 

70:  wtb  rO  ,27  ,46 , ol , int ( p2/64)  ,p2 , p3 

71:  ret 

72:  "psiz": 

73:  pi  ♦H ; p2  ♦W 

74;  wtb  r0,27,79,  int(D4*120/64)  ,p4*120,int(D3*96/64)  ,03*96 
75:  ret 
76;  "scl"; 

77:  120W/(p2-pl)-U 
78:  96H/( p4-p3)  ♦V 
79:  pl-*'X;p3^Y 
80:  ret 
81:  "xaxis"; 

82;  wtb  rO  ,27,46,95,0,5,9 

83;  if  p3=0  and  d4=0; X+p3 ;X+120M/U^p4 

34;  if  p2=0 ; P4-P3+02 

85;  wtb  r0,27,65,int((p3-X)U/64) ,int((p3-X)U) ,int( (pl-Y)V/64) 
,int((pl-Y)V) 
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86;  p3+p5;wtb  r0,43;wtb  t0,8 

87:  wtb  rO ,27 ,114 , int (p2U/64) , int ( p2U) ,0 ,0 ; wtb  r0,43,8;imp  (p 

5+p2*p5)  >=d4 

88;  ret 

89;  "yaxis"; 

90;  wtb  rO, 27, 46, 124,0 ,3,0 

91:  if  p3*0  and  p4  = 0;  Y*p3  ;y+96H/V*t54 

92;  if  p2=0 ; p4-D3*p2 

93;  wtb  r0,27,65, int ((ol-X)U/64) , int( (pl-X)U) , int( (o3-Y)V/64) 

, int( (p3-Y)V) 

94;  p3+p5;wtb  r0,43;wtb  r0,8 

95;  wtb  rO  ,27,114,0,0,  int(p2V/64)  , int(p2V)  ;wtb  rO  ,43,8  ; itno  (o 

5+p2*p5)  >=p4 

96:  ret 

97:  "space"; 

98;  if  pl<0;gto  +2 

99:  wtb  r^.,  32;jmp  2 ( (pl-l*pl)  =0) 

100:  wtb  r0,8;jinp  (pl+l*pl)»0 
101;  ret 
102:  "skip": 

103;  if  pl<0;gto  +2 

104;  wtb  r0,10;jmD  2( (pl-l*pl)  *0) 

105;  wtb  r0,27,10;jmp  (pl+l^pl)=0 
106:  ret 
107;  "form"; 

108;  wtb  r0,27,77 

109;  wtb  r0,27,84 

110;  if  pl  = 0;13. 2-*ol;ll-p2-p3 

111:  wtb  r0,27,87,int(126*pl/64)  ,120*ol 

112;  wtb  r0,27,76,int(96*p2/64)  ,96*02 

113;  wtb  r0,27,70,int(96*p3/64) ,96*03 

114;  ret 
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APPENDIX  D.  DISTRIBUTION  LIST 


NAME  AND  ADDRESS 

U,S.  Arti\y  Test  and  Evaluation  Command 
ATTN:  DRSTE-AD-M 
DRSTE-AD-I 
DRSTE-TO 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

U.S.  Arniy  Dugway  Proving  Ground 
ATTN:  STEPD-PO 
Dugway,  UT  84022 

Commander 

U.S.  Arniy  Electronic  Proving  Ground 

ATTN:  STEEP-MT-I 

Fort  Huachuca,  AZ  85613 

Commander 

U.S.  Army  Cold  Region  Test  Center 
ATTN:  STECR-TD 

STECR-TE-TD,  CPT  Cole 
APO  Seattle  98733 

Commander 

U.S.  Army  Tropic  Test  Center 

ATTN:  STETC-TD-M 

Fort  Clayton,  Canal  Zone  09827 

Commander 

U.S.  Army  Aberdeen  Proving  Ground 
ATTN : STEAP-MT-G 

Aberdeen  Proving  Ground,  MD  21005 
Commanded' 

U.S.  Army  Missile  Command 
ATTN:  DKSMI-RT 
Redstone  Arsenal , AL  35809 

Commander 

U.S.  Arn\y  White  Sands  Missile  Range 
ATTN:  STEWS- ID 
STEWS-PL 

White  Sands  Missile  Range,  NM  88002 


NO.  OF  COPIES 


3 

5 

1 


1 


1 


1 

3 


1 


1 


1 


1 

1 
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NAME  AND  ADDRESS 


NO.  OF  COPIES 


Commander 

Jefferson  Proving  Ground 

ATTN;  STEJP-TD-I  1 

Madison,  IN  47251 

Commander 

Defense  Documentation  Center 

ATTN:  Document  Service  Center  2 

Cameron  Station 
Alexandria,  VA  22314 


Commander 

U.S.  Army  Electronics  Command 

ATTN:  DRSEL-NL-P  2 

DRSEL-PD-TT  1 

DRSEL-MA  1 

DRSEL-PA  1 

DRSEL-PP-CM-C  1 

Fort  Monmouth,  NJ  07703 

Commander 

U.S.  Army  Computer  System  Command 

ATTN : CSTS-TSO-F  2 

Fort  Belvoir,  VA  22060 

President 

U.S.  Army  Aviation  Test  Board 

ATTN:  ATZQ-OT  1 

Fort  Rucker,  AL  36362 

U.S.  Army  Aircraft  Development  Test  Activity 

ATTN:  STEBG-PO  1 

Fort  Rucker,  AL  36362 
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NAME  AND  ADDRESS 


NO.  OF  COPIES 


Commander 

U.S.  Army  Yuma  Proving  Ground 

ATTN:  STEYP-MTD  1 

STEYP-MSA  3 

STEYP-MAA  1 

STYEP-MMI  25 

STEYP-MTS  1 

STEYP-MTW  1 

STEYP-MTE  I 

STYEP-MDP  1 

STEYP-MLS  1 

STEYP-TD  1 

STEYP-CO  1 

Range  Management  Office  1 

Test  Management  1 

Yuma,  AZ  85364 


t 


